Interstitial gas and density-segregation in vertically-vibrated granular media 
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We report experimental studies of the effect of interstitial gas on mass-density-segregation in a 
vertically-vibrated mixture of equal-sized bronze and glass spheres. Sufficiently strong vibration in 
the presence of interstitial gas induces vertical segregation into sharply separated bronze and glass 
layers. We find that the segregated steady state {i.e., bronze or glass layer on top) is a sensitive 
function of gas pressure and viscosity, as well as vibration frequency and amplitude. In particular, we 
identify distinct regimes of behavior that characterize the change from bronze-on-top to glass-on-top 
steady-state. 

PACS numbers: 45.70.Mg, 45.70.Qj, 64.75.+g 



Vibrated granular media exhibit characteristics both 
similar to and distinct from solids, liquids, or gases, with 
important open questions about the applicability of hy- 
drodynamic theories and the emergence of patterns and 
order 0, H & IE E| • Recent experiments have inves- 
tigated the strikingeffects of interstit ial g as on vibrated 
granular media l^lalalia lllllmlmlll . In particular, 
it was recently demonstrated [15l lid. Il7j that interstitial 
gas (in contrast to vacuum) is necessary for spatial seg- 
regation (vertical layering) by particle mass density in a 
binary granular medium of equal-sized particles that is 
vertically vibrated in a sealed cell. In this paper we re- 
port an experimental study of such density-segregation 
as a function of interstitial gas pressure and viscosity, as 
well as vibration frequency and amplitude. 

Our experimental setup was similar to that of [TH 
Hjj . The granular mixture consisted of equal-sized solid 
spheres of bronze (ACU Powder International) and soda 
lime glass (MO-SCI Corp.) with diameter d — 98 ±8 zxm, 
where 90% of both types of particles fell within this 
range. The respective mass densities of the bronze and 
glass are 8.9 g/cm 3 and 2.5 g/cm 3 . We used a volume 
mixing ratio of glass:bronze = 3:1, except where noted be- 
low. To reduce moisture, we heated each mixed granular 
sample for 3-5 minutes and used an electrically grounded 
brass pan to reduce static charge buildup. We used such 
freshly-prepared samples for each data-taking run, which 
we performed over times < 3 hours. 

We placed each sample inside a sealed, rectangular 
borosilicate glass cell (Spectrocell, Inc.) of width 46 mm, 
height 50 mm, and depth 10 mm, mounted on a vibra- 
tion platform (Fig. The height of the granular mix- 
ture at rest in the sample cell was approximately 20 mm 
for all data runs. We used an electromechanical shaker 
(Labworks, Inc. ET-126) driven by a sinusoidal wave- 
form to vibrate the cell vertically, nominally with ver- 
tical (z) position z(t) = As'm(2irft) and peak velocity 
v = 2nAf = Alo. We characterized this vibration by the 



frequency / and dimensionless acceleration T = Acu 2 / g, 
where g is the acceleration of gravity. We vibrated the 
cell over the range / = 0-200 Hz and T = 0-20. Over this 
range the particles did not interact with the top of the 
sample cell. A pair of accelerometers continuously moni- 
tored the vibration in three dimensions. We ensured that 
the dimensionless horizontal accelerations (T x and T y ) 
were small (< 0.2, except as noted below), and limited 
variation in the dimensionless vertical acceleration (r 2 ) 
to no more than 0.1 across the width of the platform. 
A gas-delivery system controlled the interstitial gas type 
and pressure inside the sample cell (P = 1-1000 Torr ± 
3 Torr), without affecting the mechanical properties of 
the cell. We performed experiments with either N2 or Ne 
gas, which have respective viscosities rj = 17.9 and 32.1 
/iPa-s at room temperature (here 77 is independent of gas 
pressure for P > 25 Torr). 

For our range of experimental conditions, the vibrated 
grains moved relative to the interstitial gas with Reynolds 
number Re « (pdv/rj) w 0.01-1, where p is the gas mass 
density. Hence viscous drag was usually larger than in- 
ertial drag. Importantly, there was no gas flow in and 
out of the sealed sample cell. Thus vertical gradients in 
the gas pressure were created across the vibrated particle 
bed inducing bulk gas flow through the granular medium 
(with gas flow velocity v g ~ VP/77 ~ P/77 as given by 
application of Darcy's Law) |7j. To avoid hysteretic ef- 
fects, we operated with the initial condition of a fully 
mixed, flat-topped sample, obtained by shaking under 
vacuum 0, E, E| ■ We then admitted the desired gas; 
observed visually and with video recording and digital 
photography (20 71m resolution) the steady-state reached 
for long-term vibration with constant Y and /, as well as 
the behavior exhibited on the approach to the steady- 
state; removed the gas; and returned the sample to its 
initial mixed configuration before another run of the ex- 
periment. 

Using this apparatus, we observed steady-state pat- 
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FIG. 1: Schematic of the apparatus, not to scale. Accelerom- 
eters are labeled by their sensitive axes (xy and xz). Dimen- 
sions are given in mm. The cell depth along the x-axis is 10 
mm (not shown) 




FIG. 2: Examples of steady-state behavior induced by ver- 
tical vibration, (a) The mixed state, which occurs under 
vacuum and was used as the initial condition for most ex- 
periments, (b) Bronze-on-top; N2 interstitial gas at P = 500 
Torr, r = 8.7, / = 90 Hz. Digital images show the boundary 
between the bronze and glass layers to be only one or two 
grain diameters thick, (c) Glass-on-top; same parameters as 
(b), except Y — 10.8. Often a glass layer also forms below the 
bronze layer, (d) Unstable pattern formation characteristic 
of very large F; N 2 interstitial gas at P = 800 Torr, V = 14, 
/ = 50 Hz. 



terns of vibration-induced density-segregation such as 
"bronze-on-top" and "glass-on-top" (see Fig. |2J) consis- 
tent with the schematic results described in [15|,|l||. We 
then studied the dependence of density-segregation and 
particle mobility on gas properties. We found that the 
heavier bronze particles are more mobile than the glass 
spheres when vertically vibrated in the presence of inter- 
stitial gas, which is consistent with particle deceleration 
due to gas-drag scaling inversely with particle mass. For 



example, we observed that vibrated bronze particles leave 
the sample surface and reach a noticeably greater height 
(about 3 times higher) than glass particles near the sur- 
face. Also, wc determined the local expansion (dilation) 
of the granular bed under vibration to be greater for re- 
gions with higher bronze concentration: images showed 
these regions to have a greater total volume than given by 
the sample's glass:bronze mixing ratio. Consistent with 
previous experiment ^(| and simulation 0], we found 
indirect evidence that the lighter glass particles are pref- 
erentially dragged along with bulk gas flow: e.g., the 
glass-on-top state never occurred when bulk gas flow due 
to gas pressure gradients was inhibited by placing the 
sample in a cell with a gas-permeable bottom |18). In 
addition, with this gas-permeable-bottom cell we found 
that overall particle mobility and bed expansion were di- 
minished. Correspondingly, T had to be increased by 
about a factor of two (by increasing the amplitude A) to 
induce bronze-on-top segregation. Further, as described 
111 all forms of density segregation — including both 
bronze-on-top and glass-on-top — were suppressed for a 
cell with a gas-permeable bottom and top. These results 
are consistent with bed expansion and particle mobility 
being enhanced by bulk gas flow through the sample. 

To characterize the competing effects of gas drag - 
i.e., greater mobility of the heavier bronze particles and 
preferential dragging of the lighter glass particles with 
bulk gas flow — we used the sealed sample cell and 
performed experiments focusing on the typically sharp 
boundary in vibration-amplitude parameter space be- 
tween the stable bronze- and glass-on-top states. In par- 
ticular, we determined the critical value of the dimcn- 
sionless acceleration (r c ) that delineates the change from 
bronze-on-top to glass-on-top steady states, as a function 
of gas pressure (P), vibration frequency (up to 200 Hz), 
and gas viscosity. Figure [21 shows examples of the mea- 
sured dependence of r c on P for typical high and low vi- 
bration frequencies (90 and 50 Hz) and gases of differing 
viscosity (Ne and N 2 ). As discussed below, we observed 
a dramatic difference in the dependence of segregation 
behavior on gas viscosity at / = 50 Hz, but little differ- 
ence at / = 90 Hz. Specifically, the glass-on-top state 
does not occur at / = 50 Hz for Ne gas at any pressure. 

Each value of r c , plotted in Fig. |3| is an average be- 
tween the highest T to yield bronze on top (Tb-hi) and 
the lowest T to yield glass-on-top (F s _; ), as determined 
by repeated sets of experimental runs. Uncertainty in 
the determination of r c comes from several sources: ir- 
reproducibility of Tb-hi and r s _; , which contributed 
an uncertainty STi w 0.1; spatial variation of Tb-hi 
and r s _; Q across the width of the platform, with asso- 
ciated uncertainty 8T2 « 0.1; and the systematic offset 
<5r3 = Tb-hi — Fg-io, which varied from w 0-0.5 for dif- 
fering experimental conditions. The total uncertainty for 
each plotted value of T c is Ar c = (STi + 6T 2 ) 1/2 + 8T3. 

^From the measurements of T c vs. P and associated 
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FIG. 3: Measured gas pressure dependence of the critical 
value of the dimensionless acceleration (F c ) delineating the 
change from the bronze-on-top to glass-on-top steady states, 
(a) / = 90 Hz vibration frequency, comparing N2 and Ne in- 
terstitial gas. (b) f — 50 Hz, N2 gas only; the glass-on-top 
state does not occur at this frequency for Ne. Four approxi- 
mate pressure regimes of characteristic r c behavior are indi- 
cated: (7) glass-on-top segregation does not occur at low gas 
pressure; (II) glass-on-top segregation occurs with the sam- 
ple approximately flat on top; {III) glass-on-top segregation 
occurs via a pressure-independent spill-over mechanism; (IV) 
glass-on-top segregation occurs via pressure-dependent spill- 
over with inertial drag nominally comparable to viscous drag. 



observations of bed dynamics we identify four charac- 
teristic regimes of steady-state behavior, as indicated in 
Fig. [3] and summarized in the following. 

Regime I (no glass-on-top steady state at low gas pres- 
sure) — The glass-on-top steady state does not occur 
at low gas pressure (below about 200 Torr for N2 and 
300 Torr for Ne), so T c is undefined. In this low pres- 
sure regime we observed minimal bed dilation and low 
mobility of the glass particles, which is consistent with 
relatively small drag from bulk gas flow driven through 
the particle bed by vibration-induced pressure gradients 
(since the gas flow velocity scales as v g ~ VP ~ P). 
Note that at very low gas pressure (< 50 Torr) we found 
that the bronze-on-top steady-state also does not occur. 

Regime II (glass-on-top occurs with minimal sample 
heaping) — At moderate gas pressures and high vibra- 
tion frequencies such as / = 90 Hz, both bronze-on-top 
and glass-on-top steady-states occur with the top of the 




FIG. 4: Time-resolved images illustrating evolution from 
an initial mixed state [like Fig. 0a)] toward a glass-on-top 
steady-state, (a) Regime II: N 2 gas, P = 250 Torr, V = 8.9, 
/ = 90 Hz. (b) Regime III: N 2 gas, P = 500 Torr, Y = 10.8, 
/ = 90 Hz. 



vibrated particle bed remaining approximately flat; i.e., 
with minimal heaping, see Fig. 0Ja) . As P is increased, 
the observed bed dilation and particle mobility increase, 
the bronze-on-top steady-state is reached more quickly 
(for r < r c ), and T c becomes larger. Similar dynamics, 
steady states, and values for T c are found for N2 and Ne 
interstitial gases. See Fig. Ota). This observed pressure- 
dependence and approximate viscosity-independence for 
r c is consistent with viscous drag due to bulk gas flow 
through the particle bed playing a dominant role [l7j . 
which implies that the drag force on the particles scales 
as F v ~ r/Vg ~ P since application of Darcy's Law gives 
v g ~ VP/rj ~ P/t) in a sealed sample cell Q- At lower 
vibration frequencies such as / = 50 Hz, Regime II be- 
havior does not occur, since for the relatively larger vi- 
bration and gas flow velocities (« Tg/cu) sample heaping 
already ensues in Regime I, i.e., at gas pressures below 
that required for glass-on-top segregation. 

Regime III (glass-on-top occurs via pressure- 
independent sample spill-over) — Above a gas pressure 
threshold, the glass-on-top state is reached via a "spill- 
over" mechanism in the presence of sample heaping, as 
illustrated in Fig. Efb). For r > T c the sample quickly 
forms a transient bronze-on-top layer, as heaping occurs 
across the width of the bed. The heaping builds until the 
bronze layer surpasses the angle of marginal stability 0] • 
Particles slide down the heap, with partial replenishment 
via convective upflow in the bronze layer. Glass particles 
then break through at the top of the heap and spill 
over, forming the glass-on-top steady-state. Particle-bed 
heaping is known to arise from vibration-induced pres- 
sure gradients and bulk gas flow, with a steepness that is 
largely pressure- and viscosity-independent for the pres- 
sures used here 0. Correspondingly, our measurements 
of T c show little dependence on P and 77 in Regime 
III at higher vibration frequencies such as / = 90 Hz. 
For lower vibration frequencies, however, we observed 
that the spill-over mechanism is inhibited relative to 
the convective replenishment of the bronze layer. In 
particular, the glass-on-top state cannot be achieved 
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for Ne gas and / < 75 Hz, nor for N2 gas and / < 40 
Hz, which is consistent with the opposite scaling of the 
viscous drag force with 77 and /: F v ~ r)V g ~ T c ij/f, 
with T c near the typical value for Regime III. 



gations of fluidized beds [2l| and may extend these stud- 
ies to segregation phenomena [22|, |23| . 

This work was supported by NSF Grant No. CTS- 
031006. 
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FIG. 5: Measured F c variation with N2 gas pressure at / = 90 
Hz, for three glass:bronze mixing ratios, 2.4:1 (squares), 3:1 
(circles), and 4:1 (triangles). As the bronze fraction grows, 
the thickness of the bronze-on-top layer increases, which sup- 
presses the mobility of the underlying glass particles and in- 
creases r c . 



Regime IV (glass-on-top occurs with pressure- 
dependent sample spill-over) At low vibration 
frequency and high gas pressure the observed difference 
in bronze and glass mobility grows with pressure, which 
inhibits the spill-over mechanism to the glass-on-top 
steady-state and yields an approximately linear de- 
pendence of T c on P. See Fig. Efb). This pressure 
dependence may result from inertial gas drag on 
the particles, which according to a simple estimate 
(Re ~ P/rjf ~ 1) becomes significant in this regime. 

Further study will be needed to establish a full, quan- 
titative understanding of the effect of interstitial gas on 
density-segregation in vertically-vibrated granular me- 
dia. For example, density-segregation can depend in de- 
tail on system parameters such as the particle mixing 
ratio. In our experiments we found that a bronze-on- 
top layer greatly suppresses the mobility of the under- 
lying glass particles; and the thicker the bronze layer 
the greater the suppression of glass mobility. One conse- 
quence of this bronze-on-top suppression is demonstrated 
in Fig. IS] introducing more (less) bronze into the system 
requires a larger (smaller) T c to bring glass to the top, 
as the glass-on-top steady-state is generally preceded by 
a transient bronze-on-top condition. In future work we 
plan to employ NMR/MRI techniques to map grain mo- 
tion as well as gas pressure and flow dynamics, by using 
NMR-detectable grains 19] and hyperpolarized noble gas 
[jfij . respectively. We are also performing NMR investi- 
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